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Abstract 



'UHE particle detection using the lunar Cherenkov technique aims to detect nanosecond pulses of Cherenkov emission 
CS| |which are produced during UHE cosmic ray and neutrino interactions in the Moon's regolith. These pulses will reach 

Earth-based telescopes dispersed, and therefore reduced in amplitude, due to their propagation through the Earth's 
D ionosphere. To maximise the received signal to noise ratio and subsequent chances of pulse detection, ionospheric 

dispersion must therefore be corrected, and since the high time resolution would require excessive data storage this 
QQ correction must be made in real time. This requires an accurate knowledge of the dispersion characteristic which is 

parameterised by the instantaneous Total Electron Content (TEC) of the ionosphere. A new method to calibrate the 

dispersive effect of the ionosphere on lunar Cherenkov pulses has been developed for the LUNASKA lunar Cherenkov 

^-p^ experiments. This method exploits radial symmetries in the distribution of the Moon's polarised emission to make 
HH Faraday rotation measurements in the visibility domain of synthesis array data (i. e. instantaneously). Faraday rotation 

measurements are then combined with geomagnetic field models to estimate the ionospheric TEC. This method of 
i-J^ ionospheric calibration is particularly attractive for the lunar Cherenkov technique as it may be used in real time to 
^p-i estimate the ionospheric TEC along a line-of-sight to the Moon and using the same radio telescope. 
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1. Introduction 

The lunar Cherenkov technique ,1] provides a promis- 
ing method of UHE neutrino detection since it utilises the 
Junar regolith as a detector; which has a far greater vol- 
ume than current ground-based detectors. This technique 
fnakes use of Earth-based radio telescopes to detect the 
coherent Cherenkov radiation emitted when a UHE neu- 
Itrino interacts in the outer layers of the Moon. It was first 
applied by Hankins, Ekers and O'SuUivan using the 64-m 
iParkes radio telescope [2|] and significant limits have been 
already been placed on the UHE neutrino flux by several 
"collaborations [l-Q- 

Electromagnetic pulses originating in the lunar surface 
will be dispersed when they arrive at Earth-based receivers 
due to propagation through the ionosphere which intro- 
duces a frequency-dependent time delay. This dispersion 
reduces the peak amplitude of a pulse, however, dedis- 
persion techniques can be used to recover the full pulse 
amplitude and consequently increase the chances of detec- 
tion. Accurate dedispersion requires an understanding of 
the ionospheric dispersion characteristic and it's effect on 
radio-wave propagation. 
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2. Effects of Ionospheric Dispersion 

The ionosphere is a weakly ionized plasma which is 
formed by ultraviolet ionizing radiation from the sun. Due 
to its relationship with the sun, the ionosphere's electron 
density experiences a strong diurnal cycle and is also de- 
pendent on the season of the year, the current phase of 
the 11-year solar cycle and the geometric latitude of ob- 
servation. The differential additive delay, caused by pulse 
dispersion, is parameterised by the ionospheric TEC (see 
Equation [ij 

At = O.OOIM X STEC X {ly-^'^ -1^-^), (1) 

where At is the duration of the dispersed pulse in seconds, 
ST EC is the Slant Total Electron Content in electrons per 
cm^ and v\o and I'hi are the receiver bandwidth band edges 
in Hz. 

Cherenkov emission produces a sub-nanosecond pulse 
and therefore detection requires gigahertz bandwidths to 
achieve the high time resolution needed to optimse the sig- 
nal to noise. Due to excessive data storage requirements, 
the only way to exploit such high data rates is to imple- 
ment real-time dedispersion and detection algorithms and 
to store potential events at full bandwidth for later pro- 
cessing. This requires an accurate knowledge of the real- 
time ionospheric TEC. 
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Ionospheric dispersion also offers some potential experi- 
mental advantages, particularly for single dish experiments 
which can not use array timing information to discrimi- 
nate against RFI. A lunar pulse will travel a one-way path 
through the ionosphere and be dispersed according to the 
current ionospheric TEC. Conversely, terrestrial RFI will 
not be dispersed at all and any Moon-bounce RFI will 
travel a two-way path through the ionosphere and be dis- 
persed according to twice the current ionospheric TEC. 
Therefore performing real-time ionospheric dedispersion 
will optimise detection for lunar pulses and provide dis- 
crimination against RFI. Dispersion may also be seen to 
offer an increase in dynamic range. If triggering is per- 
formed on a dedispersed data stream while the raw data 
is buffered, any pulse clipping that occurs in the trigger- 
ing hardware can be recovered during offline processing by 
reconstructing the pulse from the raw, undispersed data. 

3. Dedispersion Hardware 

Pulse dispersion can be corrected using matched filtering 
techniques implemented in either analog or digital technol- 
ogy. Early LUNASKA experiments made use of the Aus- 
tralia Telescope Compact Array (ATCA) which consists 
of six 22-m dish antennas. Three antennas were fitted 
with custom-designed hardware for the neutrino detection 
experiments and pulse dedispersion was achieved through 
the use of innovative new analog dedispersion filters that 
employ a method of planar microwave filter design based 
on inverse scattering [7] . 

As the microwave dedispersion filters have a fixed dedis- 
persion characteristic, an estimate had to be obtained for 
the TEC which would minimise errors introduced by tem- 
poral ionospheric fiuctuations. The ATCA detection ex- 
periments were performed in 2007 and 2008 during solar 
minimum and therefore relatively stable ionospheric condi- 
tions. Initial observations were during the nights of May 5, 
6 and 7, 2007 and these dates were chosen to ensure that 
the Moon was at high elevation (particularly during the 
night-time hours of ionospheric stability) and positioned 
such that the ATCA would be sensitive to UHE particles 
from the galactic center. The filter design was based on 
predictions made using dual-frequency GPS data and as- 
sumed a differential delay of 5 ns across the 1.2-1.8 GHz 
bandwidth. Data available post experiment revealed that 
the average differential delay for these nights was actually 
4.39 ns, with a standard deviation of 1.52 ns. 

The ionosphere experiences both temporal and spatial 
fluctuations in TEC and therefore some signal loss is ex- 
pected with a fixed dedispersion filter. A promising dig- 
ital solution to overcome these losses lies in the use of 
high speed Field Programmable Gate Arrays (FPGAs). 
An FPGA implementation allows the dedispersion charac- 
teristic to be tuned in real time to reflect temporal changes 
in the ionospheric TEC. A fully coherent or predetection 
dedispsersion method was pioneered by Hankins and Rick- 
ett which completely eliminates the effect of disper- 



sive smearing. This is achieved by passing the predetected 
signal through an inverse ionosphere filter which can be 
implemented in either the time domain, as an FIR filter, 
or in the frequency domain. 

In 2009, the LUNASKA collaboration started a series 
of UHE neutrino detection experiments using the 64-m 
Parkes radio telescope. For these experiments, dedisper- 
sion was achieved via a suite of FIR filters implemented on 
a Vertex 4 FPGA. As GPS TEC estimates are currently 
not available in real-time, near real-time TEC measure- 
ments were derived from foF2 ionosonde measurements, 
lonosondes probe the peak transmission frequency (fo) 
through the F2-layer of the ionosphere which is related to 
the ionospheric TEC squared. A comparison to GPS data 
available post-experiment revealed that the foF2-derived 
TEC data consistently underestimated the GPS TEC mea- 
surements. This is attributed to the ground-based lonoson- 
des probing mainly the lower ionospheric layers and not 
properly measuring TEC contribution from the plasmas- 
phere [10|. 

4. Monitoring the Ionospheric TEC 

Coherent pulse dedispersion requires an accurate knowl- 
edge of the ionospheric dispersion characteristic which is 
parameterised by the instantaneous ionospheric TEC. 

TEC Measurements can be derived from dual-frequency 
GPS signals and are available online from NASA's CDDIS 
11| . however, these values are not available in real time. 
The CDDIS TEC data is sampled at two-hour intervals 
and is currently published after at least a few days delay. 
Estimates derived from foF2 ionosonde measurements are 
available hourly from the Australian Ionospheric Predic- 
tion Service 12|. However, as discussed, there are known 
inaccuracies in the derivation of the foF2-based TEC esti- 
mates. 

Both of these products are published as vertical TEC 
(VTEC) maps which must be converted to Slant TEC 
(STEC) estimates to obtain the true total electron con- 
tent through the slant angle line-of-sight to the Moon. To 
perform this conversion, the ionosphere can be modeled as 
a Single Layer Model (SLM) [l^ which assumes all free 
electrons are concentrated in an infinitesimally thin shell 
and removes the need for integration through the iono- 
sphere. Slant and vertical TEC are related via 

STEC ^ F{z)VTEC. (2) 
where F{z) is a slant angle factor defined as 
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Re is the radius of the Earth, z is the zenith angle to the 
source and H is the height of the ideahsed layer above 
the Earth's surface (see Figure [IJ . The CDDIS also use 
an SLM ionosphere for GPS interpolation algorithms and 
assume a mean ionospheric height of 450 km. 
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Figure 1: Parameters of the Ionospheric Single Layer Model. 



5. A New Method of Ionospheric Calibration 

As the solar cycle enters a more active phase accurate 
pulse dedispersion is becoming a more important experi- 
mental concern for the lunar Cherenkov technique. This 
requires methods of obtaining more accurate measure- 
ments of the ionospheric TEC. 

A new technique has been formulated to obtain TEC 
measurements that are both instantaneous and line-of- 
sight to the Moon. The ionospheric TEC can be deduced 
from Faraday rotation measurements of a polarised source 
combined with geomagnetic field models, which are more 
stable than ionospheric models (the CCDIS [ll| states 
that ionospheric TEC values are accurate to '^20% while 
the IGRF magnetic field values are accurate to bet- 
ter than 0.01%). Lunar thermal emission can be used as 
the polarised source since Brewster angle effects produce 
a nettpolarisation excess in the emission from the lunar 
limb [15| . This provides a method for calibrating the iono- 
sphere directly line-of-sight to the Moon and makes the 
lunar Cherenkov technique extremely attractive for UHE 
cosmic ray and neutrino astronomy as it allows the char- 
acteristic dispersion to be used as a powerful discriminant 
against terrestrial RFI whilst removing the need to search 
in dispersion-space. 

The unique constraints of an UHE neutrino detection 
experiment using the lunar Cherenkov technique conflict 
with traditional methods of planetary synthesis imaging 
and polarimetry which requires a complete set of spacings 
and enough observing time for earth rotation. Therefore to 
apply this method of ionospheric calibration to the ATCA 
detection experiments, innovations in the analysis of lunar 
polarisation observations were required. In particular, a 
method of obtaining lunar Faraday rotation estimates in 
the visibility domain (i. e. without Fourier inversion to 
the image plane) had to be developed. 



Working in the visibility domain removes both the imag- 
ing requirement of a compact array configuration, which 
would increase the amount of correlated lunar noise be- 
tween receivers, and also removes the need for earth rota- 
tion allowing measurements to obtained in real time. This 
technique makes use of the angular symmetry in planetary 
polarisation distribution. The intrinsic thermal radiation 
of a planetary object appears increasingly polarised toward 
the limb, when viewed from a distant point such as Earth 



151 . Il6j . The polarised emission is radially aligned and is 
due to the changing angle of the planetary surface toward 
the limb combined with Brewster angle effects. The an- 
gular symmetry of this distribution can be exploited by 
an interferometer so that an angular spatial filtering tech- 
nique may be used to obtain real-time position angle mea- 
surements directly in the visibility domain. The measured 
position angles are uniquely related to the corresponding 
uv angle at the time of the observation. Comparison with 
the expected radial position angles, given the current uv 
angle of the observation, gives an estimate of the Faraday 
rotation induced on the Moon's polarised emission. Fara- 
day rotation estimates can be combined with geomagnetic 
field models to determine the associated ionospheric TEC 
and subsequently provide a method of calibrating the cur- 
rent atmospheric effects on potential Cherenkov pulses. 

Observations of the Moon were taken using the 22-m 
telescopes of the Australia Telescope Compact Array with 
a center frequency of 1384 MHz. At this frequency the 
Moon is in the near field of the array, however, investiga- 
tion of the Fresnel factor in polar coordinates showed that 
it has no dependence on the spatial parameter, which de- 
termines the polarisation distribution of a planetary body. 

Using the angular spatial filtering technique, position 
angle estimates were calculated directly in the visibility 
domain of the lunar observational data. The Faraday rota- 
tion estimates were obtained by comparing these angles to 
the instantaneous uv angle and the resultant Faraday ro- 
tation estimates were averaged over small time increments 
to smooth out noise-like fiuctuations. Since the polarised 
lunar emission received on each baseline varied in intensity 
over time, there were nulls during which the obtained po- 
sition angle information was not meaningful. A threshold 
was applied to remove position angle measurements taken 
during these periods of low polarised intensity and base- 
line averaging was considered necessary as the results on 
each baseline were slightly different and each affected dif- 
ferently by intensity nulls. The Faraday rotation estimates 
were converted to estimates of ionospheric TEC via 



n = 2.36 X lO^iy- 



N{s)B{s)cos{e)ds, (5) 



path 



where fl is the rotation angle in radians, / is the signal 
frequency in Hz, N is the electron density in m'^, B is the 
geomagnetic field strength in T, 6* is the angle between the 
direction of propagation and the magnetic field vector and 
ds is a path element in m. 
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Figure 2: Lunar Faraday rotation estimates converted to (left) iono- 
spheric TEC values and (right) the differential delay across 1.2-1.8 
GHz. 



To evaluate the effectiveness of this new ionospheric cah- 
bration technique, the TEC resuhs were compared against 
ionospheric TEC estimates derived from dual-frequency 
GPS data (Figure ^ . Slant angle factors were used to 
convert the GPS VTEC estimates to STEC toward the 
Moon for comparison with the ATCA data. Both data 
sets exhibited a similar general trend of symmetry around 
the Moon's transit. However, the ATCA data often under- 
estimated the GPS data, particularly around 14:30-17:00 
UT where the STEC estimates may have been influenced 
by bad data from the shorter baselines or due to TEC 
contribution from the plasn iasp here which is not in the 
presence of a magnetic field [lO|. These observations were 
taken when the TEC was very low and therefore the rela- 
tive error in the TEC estimates is large. 

6. Conclusions 

As the sun enters a more active phase, accurate iono- 
spheric pulse dispersion is becoming a more important ex- 
perimental concern for UHE neutrino detection using the 
lunar Cherenkov technique. Hardware dedispersion op- 
tions rely on the accuracy of real-time ionospheric TEC 
measurements and, while there are a few options avail- 
able for obtaining these measurements, they are not cur- 
rently available in real time nor directly line-of-sight to the 
Moon. A new ionospheric calibration technique has been 
developed. This technique uses Faraday rotation mea- 
surements of the polarised thermal radio emission from 
the lunar limb combined with geomagnetic field models to 
obtain estimates of the ionospheric TEC which are both 
instantaneous and line-of-sight to the Moon. STEC esti- 
mates obtained using this technique have been compared 
to dual-frequency GPS data. Both data sets exhibited sim- 
ilar features which can be attributed to ionospheric events. 



however, more observations are required to investigate this 
technique further. 
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